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REGRESSION VERSION OF THE MATSUMOTO-YOR TYPE 
CHARACTERIZATION OF THE GAMMA AND KUMMER DISTRIBUTIONS 

JACEK WESOLOWSKI 


Abstract. In this paper we study a Matsumoto-Yor type property for the gamma and Rummer inde¬ 
pendent variables discovered in Koudou and Vallois (2012). We prove that constancy of regressions of 
U = (1 + [X + Y) -1 )/(l + A' -1 ) given V = X + Y and of t/ _1 given V, where X and Y are indepen¬ 
dent and positive random variables, characterizes the gamma and Rummer distributions. This result 
completes characterizations by independence of U and V obtained, under smoothness assumptions for 
densities, in Koudou and Vallois (2011, 2012). Since we work with differential equations for the Laplace 
transforms, no density assumptions are needed. 


1. Introduction 

Let X and Y be independent random variables. There are several well known settings in which 
U = if>(X,Y) and V = X + Y are also independent. Related characterizations of distributions of X and 
Y by properties of independence of X and Y and independence of U and V have been also studied. The 
most prominent seem to be: 

• Bernstein (1941) characterization of the normal law by independence oiU = X — Y and V, 

• Lukacs (1956) characterization of the gamma law by independence of U = X/Y and V. 

In the end of 1990’s a new independence phenomenon of this kind, called Matsumoto-Yor property, see 
e.g. Stirzaker (2005), p. 43, was discovered. It says that for X with a GIG (generalized inverse Gaussian) 
law and independent Y with a gamma law (both distributions with suitably adjusted parameters), random 
variables U = 1/X — 1/(X + Y) and V are independent. This elementary property was identified while 
the authors analyzed structure of functionals of Brownian motion - see Matsumoto and Yor (2001, 2003). 
A related characterization of the GIG and gamma laws by independence of X and Y and of U and V 
was obtained in Letac and Wesolowski (2000). Both these results: the Matsumoto-Yor property and 
the characterization were generalized in several directions. Matrix variate analogues were studied e.g. 
in Letac and Wesolowski (2000), Wesolowski (2002) and Massam and Wesolowski (2006) - the last one 
including a relation with conditional structure of Wishart matrices. Recently it has been extended to 
symmetric cones setting in Kolodziejek (2014). Multivariate versions related to specific transformations 
governed by a tree were considered in Barndorff-Nielsen and Koudou (1998), Massam and Wesolowski 
(2004) and Koudou (2006) . Further connections with (exponential) Brownian motion were investigated 
in Wesolowski and Witkowski (2007) and Matsumoto et al. (2009). There are also regression versions 
of Matsumoto-Yor typ characterizations, as given in Seshadri and Wesolowski (2001), Wesolowski (2002) 
and Chou and Huang (2004). A survey of these results together with other characterizations of the GIG 
law can be found in a recent paper Koudou and Ley (2014). 

In 2009 Koudou and Vallois tried to generalize Matsumoto-Yor property by a search of distributions 
of independent X and Y and functions / such that V = f(X + Y) and U = f{X) — f(X + Y ) are 
independent. Their research lead to a discovery of another pair U = ijj(X,Y) and V = X + Y with 
independence property: Assume that X and Y are independent random variables, X has the Kummer 
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distribution K(a, b, c) with the density 

fx(x) oc I( 0 <oo )(x), a,b,c> 0, 

and Y has the gamma distribution G( 6 , c) with the density 

fy{y) oc y b ~ 1 e-- cy I(o,oo) {y)- 

Then, see Koudou and Vallois (2012), random variables 

U = 1 +X t^ and V = X + Y 


( 1 ) 


i+ 


x 


are independent, U has the beta first kind distribution B/(a, b) with the density 


fu(u)cxu a l {l-u) b 1 / ( 0 i 1 ) ( m ) 

and V has the Kummer distribution, K(a + 6 , — 6 , c). (Note that the Kummer distribution K(a,f3, 7 ) is 
well-defined iff a, 7 > 0 and /3 G M.) 

It is an interesting question if a theory, similar to the one for the original Matsumoto-Yor property 
described in the literature recalled above, can be developed for this new independence property. There 
have already been some successful efforts in this direction. The property was extended to matrix variate 
distributions in Koudou (2012). It is also known, see Koudou and Vallois (2011, 2012), that, under 
appropriate smoothness assumptions on densities, a characterization counterpart of the property holds: 
if X and Y are independent positive random variables, and U and V, given by m, are also independent 
then X ~ K(a, 6 , c) and Y ~ G(b, c) for some positive constants a, 6 , c. These smoothness restrictions 
require that the densities of X and Y are strictly positive on (0, 00 ) and either twice differentiable or their 
logarithms are locally integrable. Letac (2009) conjectured that such a characterization is possibly true 
with no assumptions on densities. In this note we contribute further to this development following the 
characterization path. Our aim is to show a characterization of the gamma and Kummer laws without 
referring to densities at all. Actually, we will consider constancy of regressions condition which, up to 
necessary moment assumption, is weaker than independence. 


2. Regression characterization 


Our main result is a characterization of the Kummer and gamma laws by constancy of regressions of U 
and I / -1 given V in the setting described in ©■ Since U G (0,1) P-a.s. EU < 00, and one can consider 
conditional moment E(t/|V) without any additional restrictions. This is not the case of E(t/ _ 1 |V) since, 
a priori, the moment EC / -1 may not be finite. Since 


u - 1 


1+-Y X+Y < 1 , J_ 
X 1 +X+Y — 1 ' X 


we have E/7 -1 < 1 +EX -1 . So, under the assumption EX -1 < 00 the conditional moment E(I/ _ 1 |K) 
is well defined. 

Now we are ready to state the main result of this note. 


Theorem 2.1. Let X and Y be independent positive non-degenerate random variables and E X 1 < 00 . 
Define U and V through ©• If 

( 2 ) E(J7|V)=a and E(t/- 1 |V) = fi 


for real constants a and then there exists a constant c > 0 such that 


X ~K 


(l + ^, (1 a) a ( P 1 } ,c) and 


Y ~ G 


/ (l—a)(/3 — 1) \ 

V «/3 ’ C J • 
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Proof. First, rewrite the equations 0 as 


E 


(-2L- 
V i+x 


X + Y \ = a - Y+Y 

^ ^ 1 1 +X+Y 


and E(±±A|X + y) 


Equivalently, we have 

( 3 ) 

and 


E (irr|* + y ) 


X+Y =l-a 


1+X+Y 


a i+x+y 
y X+Y 


( 4 ) 

The equation © implies 

( 5 ) 


E (Tl^ + ^)=/3-l + Tf 


E 


e »( l + X + V) 

T+X 


(l-a)Ee s(1+x+y > + cdE 


e «(l +X+Y) 

l+X+Y 


at least for s < 0 . 

Similarly, from 0 we get the equation 
( 6 ) E eS( '^ +ri = ((3 - l)Ee s(x+y) + /3E e ^T ) , s < 0. 

Differentiating 0 with respect to s (it is possible at least for s < 0) we obtain 


EeS (i+Ar+r) +E y e s{x+Y+i) = (i - a )E (1 + X + F)e s(1+x+y) +aEe s(x+r+1) . 


After dividing by e s both sides of the above equation and canceling the term Ee s ^ A+1 ^ we arrive at 
E E Ye sY = (1 - a) (E Xe sX E e sY + E e sA E Ye sV ) . 

This equation can be written as 

(7) e~ s K M' = (1 — a)(LM)', 


where 

Jf(s)=E^, L(s) = Ee sX and M{s) = Ee sY . 

On the other hand differentiating 0 we get 

Ee s(x+y) + E y e s(X+Y) = (/3 - 1 )E (X + Y)e s(x+Y) + /3Ee s(A ' +Y) . 

Consequently, 

E ^ E Ye sY = (/3- 1) (E Xe sA E e sY + E e sA E Ye sY + E e sA E e sY ) . 

Therefore 


( 8 ) GM' = (/3-l)((LM)'+ LM), 

where 

G(s) = E^. 

By deriving the formula for (. LM)' from 0 and 0 get 

(9) ae~ s KM' = bGM' - LM, 
with a = (1 — a ) -1 and b = (0 — l) -1 . Differentiate 0 to get 

- ae~ s KM' + ae~ s K'M' + ae~ s KM" = bG'M' + bGM" - (LM)'. 

Note that G' = L = e~ s K'. Therefore the above equation together with 0 and 0, after multiplying 
both sides by M' implies 

— (LM)'M' + aLM' 2 + (LM)'M" = bLM' 2 + ((LM)' + LM)M" - (LM)'M' 
which after cancelations (which are allowed in a left neighborhood of zero) gives 

MM" = (a- b)M' 2 . 
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Note that 


a — h— _J |—L_ = 2 ~ q ~P =1 + 

“ U — 1 — ol t 1—0 — (1—q)(1—/3) — 1 + 


a/3-1 _1 4- 1 

(1 —a)(/3 —1) — 1 + V 


and due to obvious inequalities: a < 1, /3 > 1 and a/3 > 1, we conclude that p > 0. Consequently, Y has 
a gamma distribution G{p, c), where c is a positive constant. 

Now we differentiate equation © for s < 0 getting 

bG'M' + bGM" = (. LM)" + {LM)'. 

Multiplying both sides by M' and using again © we arrive at 

bLM' 2 + {{LM)' + LM)M" = {LM)" M' + {LM)'M' 

which yields 


T H M I T I 

L M + L 


= o. 


(M'\ , M' M" t M" (M' \" 

\M ) + M M h M l 1 b > \ M ) 

After inserting known values for M, M' and M" the above equation transforms into 

(c — s)L"{s) + (1 — p + c — s)L'{s) — (1 + bp)L{s) = 0, s < 0. 

Change the variable t c — s and define N{t) = L{c — t). It follows that 

tN"{t) + (1 — p — t) N'{t) — (1 + bp) N{t) = 0 t > c. 

We read two linearly independent solutions of this equation from Abramovitz and Stegun (1965), Ch. 
13. One of these solutions is the generalized hypergeometric function 

N{t) = M{ 1 + bp , 1 -p,t) =iFi(l + bp , 1 -p,t) 

which is of the order e t t^ 1+p ^ b for t -> oo, see (13.1.4) in Abramovitz and Stegun (1965), and thus tends 
to infinity with s —» — oo and thus t = c — s —> oo. In the case we consider this is impossible since the 
Laplace transform of negative argument s of positive probability measure has to be bounded. The second 
solution 


N{t) = U{\ + bp , 1 — p,t) = C f 

Jo 


e ~ tx __ dx 


yields 


L{s) = C f 

Jo 


(14z) p ( b+1 )+i ( 


dx, 


which is a Laplace transform of the Rummer K {bp + 1 ,p, c) distribution. 


□ 


Remark 2.1. Recall that U ~ Beta/(1 + bp,p) and thus EX 1 < oo, as expected. Moreover, since p > 1 
then also E (1 — /7) -1 < oo. 

Remark 2.2. R still not clear if independence ofU and V for independent, positive and non-degenerate 
X and Y without any additional assumptions characterizes the gamma and Rummer laws. Theorem 1 
answers the question under additional restriction that EU -1 < oo. 

Remark 2.3. Since U as defined in CD is (0,1) valued random variable, without any additional moment 
assumptions we can write regressions conditions of the form 

E ((1 — U) k \ V) = atk 

for some positive k’s and 04 - ’s. Obviously, such conditions are weaker than independence. A little of 
algebra allows to see that the above condition is equivalent to 


E 


((^jx + r) =a k {X + Y) k . 


However, we failed to prove characterization assuming the above conditions for, say, k = 1,2. 
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